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A modified smoothed particle hydrodynamic (MSPH) computational technique was utilized to simulate 
molten particle motion and infiltration speed on multi-scale analysis levels. The radial velocity and velocity 
gradient of molten alumina, iron infiltration in the TiC product and solidification rate, were predicted 
during centrifugal self-propagating high-temperature synthesis (SHS) simulation, which assisted the 
coating process by MSPH. The effects of particle size and temperature on infiltration and solidification of 
iron and alumina were mainly investigated. The obtained results were validated with experimental 
microstructure evidence. The simulation model successfully describes the magnitude of iron and alumina 
diffusion in a centrifugal thermite SHS and Ti + C hybrid reaction under centrifugal acceleration. 

Self-propagating high-temperature synthesis (SHS), or combustion synthesis (CS), is an advanced method of 
producing high-temperature materials, such as refractory or intermetallic materials and cermets. The CS 
process is characterized by extreme heating rates, short reaction times, and expelling impurities during 
reaction 1 ' 2 . Ceramic-lined steel pipes are a sort of high-performance, heat-resistant, anti- corrosion and wear- 
resistant pipe that can be produced by a centrifugally- assisted thermite (CT) process 3 . The ceramic-lined steel 
pipe produced by SHS technique has occupied the market rapidly. Not only are the quality superior and the 
capabilities outstanding, but the price ratio is higher than of other wearable, anti-abrasion and anti-heat pipes. 
Moreover, they weigh less compared to other wearable, cast alloy steel pipes 4 ' 5 . 

In a research by Gao, et al. 6 , in China, W-C-Fe cermet-lined steel pipes were produced through a SHS- 
centrifugal process. Finite element analysis of residual thermal stress was performed on ceramic-lined composite 
pipes prepared by centrifugal SHS. The finite element method of simulating the temperature and residual thermal 
stress distribution of the ceramic-lined composite pipe made by centrifugal SHS was discussed. Wang and Yang 7 
also employed various SHS layer thicknesses to fabricate centrifugal- SHS composite pipes with the aim of 
reducing residual stress. Therefore, the reliability of the ceramic-lined composite pipes should improve. 

Local reinforcement of functionally graded coating produced using a centrifugal-assisted thermite process has 
recently come to the attention of researchers 8 " 10 . The thermite reaction (A1 2 0 3 and Fe) infiltrated the TiC pellets 
and created a strong, titanium aluminide intermetallic layer. There is a compelling need to simulate the mech- 
anism of locally reinforcing lined ceramic tubes. Several simulation techniques are applied in SHS, among which 
finite difference numerical method 11 , numerical modeling of field- activated combustion synthesis process 12 , 
modeling the solidification of functionally graded materials by centrifugal casting 1013 , and analytical 
modeling 1415 . 

Initially developed inl977 1617 , Smoothed Particle Hydrodynamics (SPH) is a meshless Lagrangian particle 
method for obtaining approximate numerical solutions of fluid dynamics equations. The SPH method was 
originally meant for astrophysical applications. It gradually extended to problems concerning incompressible 
fluids using either weakly compressible fluid models, or algorithms designed to solve full incompressible equa- 
tions 18 . In this method, the fluid is replaced by a set of particles 19 ' 20 . At present, SPH has a wide range of 
applications in numerous fields, such as heat conduction, gas explosions, micromachining, micro -forming, 
and granular flows that have demonstrated the ability to simulate highly non-linear free-surface flows including 
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Figure 1 | FESEM microstructure of functionally graded TiC-Fe-Al203 
coating; (a) Region next to the inner most layers relative to the pipe's axis: 
zone (I) Alumina-rich layer coating on the Ti + C pellet, (II) Ti-Fe-Al-rich 
zone intermetallic on the Ti + C pellet and the coated layer, (III) Fe-rich 
zone in the Ti + C pellet; (b) a higher magnification of zone (II). 



wave overturning, jets, and the formation of spray and droplets 21,22 . 
SPH has been vastly developed for different applications like under- 
water explosion 20 , wave propagation 19 , magnetohydrodynamics 23 , 
bulk deformation 24 , and particle motion 18 . SPH has also been imple- 
mented to solve complex splashing free- surface flows and the dif- 
ferential motion of multiple solid-casting processes owing to its 
mesh-free nature. It can additionally handle boundaries and apply 
fragmenting free-surface flow of solid particles 25 . 

In the SPH technique, the spatial gradients of speed and deposition 
rate at the boundaries are difficult to calculate due to insufficient 
data. This method also employs ghost particles at the boundary, 
which are not real or adequate for solving approximation problems 
at the boundary. Therefore, to mitigate such drawbacks, modified 
smoothed particle hydrodynamics (MSPH) is adopted based on the 
Taylor expansion as opposed to ghost particles 26 . 

Mathematically, the range of problems that can be solved by 
MSPH is much broader than by mesh-based methods, since approxi- 
mation is not based on elements whose distortion may degrade cal- 
culation accuracy, something helpful in both the fluid and solid 
domains 27 . However, to the best of the authors' knowledge, there 
are no reports on meshless particle methods such as MSPH to simu- 
late SHS particle behavior, e.g., velocity, viscosity, and displacement. 

Therefore, the purpose of this work is to implement MSPH to 
predict the radial velocity and velocity gradient of molten alumina, 
iron infiltration in the TiC product and solidification rate during the 
centrifugal SHS-assisted functional coating process. The effects of 
particle size and temperature on the infiltration and solidification 
of iron and alumina are mainly investigated. The obtained results are 
validated with experimental evidence of microstructure. 

Results 

Temperature data obtained from the solution of the heat balance 
equation was employed to calculate the fluids' (iron and alumina) 
viscosity, particle speed in the radial direction, gradient of speed, and 
particles' position. A microstructural analysis is also presented in this 
section. In order to verify the MSPH code, the exact solution of 
velocity equation (7) and values obtained from the MSPH solution 
are plotted accordingly. In addition, the measured temperatures are 
compared with the MSPH solution results from equation (6). 

Microstructural analysis. A typical area of functionally graded 
TiC-Fe-Alumina functionally graded coating, processed using 
centrifugally assisted combustion synthesis method, was taken for 
microstructural observations. A field emission scanning electron 
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microscopy (FESEM) micrograph of the typical area is shown in 
Figure 1. Figure 1 (a) illustrates a region next to the pellet's inner 
surface relative to the rotation axis, while Figure 1 (b) depicts a higher 
magnification of zone (III) from Figure 1(a). Corresponding energy 
dispersive spectroscopy (EDS) elemental analysis of Figure 1 (a) for 
zones (I), (II), and (III) is listed in Table 1. Zones (I), (II), and (III) in 
Figure 1 (a) are marked relative to the tube's axis of rotation. 

In accordance with the EDS and FESEM results, the thickness of 
the alumina-rich layer comprising zones (I) and (II) is roughly 
110 |im. Moreover, zone (I) is a composite- coated layer with an 
alumina-rich constituent. However, based on micrograph measure- 
ments, alumina diffusion occurs around zone (II), whose thickness is 
around 30 |im. On the other hand, iron particles diffuse deeper into 
the Ti + C pellet compared with alumina. 

Particle viscosity versus time. The time-viscosity plot of a centri- 
fugal thermite reaction is illustrated in Figure 2 along with the corre- 
sponding temperature cooling reaction. The time-viscosity plot is the 
MSPH solution of the exact calculations using Arrhenius' equation 
adopted in eq. (8). 

Particle velocity and velocity gradient. The exact particle solution is 
derived based on fundamentals of fluid mechanics, which provides a 
closed-form solution for the velocity and position of alumina and 
iron as given by equations (7) and (11). The exact solution of particle 
velocity and the corresponding MSPH simulation in three different 
particle size cases (dp, dp2, and dp3) are plotted in Figure 3. The 
simulation parameters are listed in Table 2. 

Figure 4 shows the radial velocity of molten alumina particles for 
three different particle sizes as predicted by MSPH and the exact 
solution. Figure 5 illustrates the particle velocity gradient in the radial 
direction as estimated by equation (11) based on the MSPH model. 

Change in particle position. The instantaneous change in radial 
particle position versus time during the deposition process is calcu- 
lated for iron and alumina particles and the results are given in 
Figure 6. Evidently, there is an indirect correlation with viscosity 
change over time. 

Discussion 

Once the thermite reaction occurred under the influence of centri- 
fugal acceleration, the molten Al 2 0 3 -Fe thermite products infiltrated 
the Ti + C pellet. The infrared pyrometer recorded the reaction 
temperature released up to 2830°C (Figure 2). The thermite molten 
products were pushed to the inner surface of the Ti + C pellet under 
the applied centrifugal force. The molten semi-products were for- 
cibly trying to penetrate the porous Ti + C media to form an in-situ 
TiC-Fe-Al 2 0 3 composite. According to the temperature observation, 
the reaction lasted around 2.5 seconds from ignition until Fe solidi- 
fication, at which point the reaction and infiltration presumably 
stopped. The thermite reaction products (Fe and A1 2 0 3 ) deposited 
according to their densities 28 at the graphite mold's innermost layer, 
after which they diffused into the TiC zone. The reactions created a 
porous TiC product. A more in-depth discussion on the current 
method is given elsewhere 8 . 

In Figure 1 (a), zone (I) is alumina-rich, and zones (II) and (III) are 
Fe-rich. Zone (I) is located at the innermost layer relative to the pipe's 
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axis. Alumina mainly formed a coating layer on the TiC pellet and 
slightly diffused into it, whereas the iron phase mostly diffused inside 
the TiC pellet. The composition of the phases gradually changed with 
the specimen's volume. Figure 1 (b) represents the magnified micro- 
graph corresponding to zone (III). It is clear how the molten iron 
diffused inside the porous area of the compacted Ti + C pellet and 
formed an interesting composite. Iron particles exist in all layers, 
with a declining rate from the inner to the outermost layer relative 
to the axis of rotation. Centrifugal acceleration significantly affected 
both metallurgical alloying and mechanical interlocking between 
different specimen layers during product formation. It is apparent 



that the alumina phase did not significantly diffuse in zone (III), 
whereas the iron phase diffused into the TiC pellet. 

As revealed in Figure 2, the viscosity of molten alumina particles 
exhibits higher values and abrupt changes at 0.8 and 2.4 s after the 
reaction started, compared to iron particles. This is because the melt- 
ing point of alumina is above 2200°C and it starts to solidify at 0.7 s, 
while Fe is in a molten state seeing that its melting point is 1535°C. 
This explains why Fe was able to penetrate deeper into TiC, as well as 
according to the FESEM microstructure observation of Figure 1, 
zone (iii). The viscosity analysis serves to calculate the particle velo- 
city and location at any given time. 
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Figure 3 | Iron particle radial velocity versus time obtained from the analytical model and MSPH predictions for different particle sizes. 
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With respect to Figure 3, in the early process stages and shortly 
after the reaction began, the iron particles' radial velocity slowed 
down rapidly as viscosity increased. It completely stopped after three 
seconds of reaction process. Physically, the particle migration rate 
(particle velocity) becomes zero when the viscosity sharply increases 
during solidification or when the particles reach a boundary that 
halts further movement. Under the centrifugal field where the part- 
icles were larger, the initial speed would be higher, compared with the 
smaller particles. The particle velocity of MSPH prediction is in close 
correlation with the exact solution results from Figure 3. This vali- 
dates the simulation model and permits for the application of MSPH 
in the SHS process. As seen in Figure 4, the alumina particle velocity 
seems to be approximately half the velocity of iron particles for 
corresponding particle sizes. This may be attributable to the density 
of iron, which is about twice that of alumina. 

Equation (1 1) is simpler than deriving and using the exact express- 
ion for velocity gradient calculations. Obviously, the velocity gra- 
dient can describe continuous particle deceleration for particles of 
different sizes. The exponential deceleration of iron particles was 
predicted because iron was in a molten state for longer than alumina, 
as observed from the temperature-time curve of the thermite reac- 
tion. With increasing viscosity (Figure 2), particle motion would 
naturally be affected. For instance, in Figure 5, as the particle size 
increases, the velocity gradient slope is greater and the particles 
decelerate faster. However, this does not drastically happen for alu- 
mina, whereby the particle speed gradient curve tends to become 
zero almost at the same time regardless of particle size. Never- 
theless, the gradient of velocity in the radial direction for iron sig- 
nificantly correlates to particle size. 

When the velocity gradient becomes zero at this stage, there might 
not be molten phase diffusion; but the solid-state diffusion phenom- 
enon is probably still mainly happening 29,30 for iron particles. Solid- 
state diffusion is widely accepted with respect to combustion 
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Figure 4 | Alumina particle radial velocity (for three particle sizes) and 
the corresponding MSPH simulation. 



synthesis processing. When the material has less motion in the liquid 
phase and owing to the in -situ high temperature gradient, the part- 
icles diffuse into each other and form intermetallics, alloys or cer- 
mets. If this phenomenon combines with centrifugal force, there 
would be extra external force acting beyond the viscosity increment, 
which poses resistance to the particles' motion. Finally, there would 
be a greater amount of the obtained product than expected. Velocity 
and velocity gradient results are in agreement with the experimental 
evidence of the centrifugal thermite and molten iron infiltration in 
the TiC pellet discussed earlier 8 . Moreover, the simulation results 
help infer the reason why alumina particles are observed at the TiC 
pellet boundary and why iron particles infiltrated the TiC porous 
medium (Figure 1). 

The indirect correlation with viscosity change over time is mani- 
fest in Figure 6. The particle position becomes zero exactly when the 
velocity gradient becomes zero. The particle motion completely stops 
after the solidification process takes place for 3 seconds. Sudden 
changes in molten particle position are obvious when the viscosity 
decreases, and the particles have no relative motion and displace- 
ment after 0.7 s from the SHS reaction. However, the particles may 
experience solid diffusion. Consequently, the particle diffusion speed 
tends toward zero as time increases. 

Again, it is understood that alumina particles are not moving 
toward the TiC pellet as the iron particles are (Figure 6). The alumina 
particles are mostly deposited at the internal boundary of the TiC 
pellet (Figure 1), so-called the top surface, due to the rapid increase in 
relative viscosity (Figure 2). Meanwhile, the viscosity of alumina 
increases in a fraction of a second, whereas the iron viscosity does 
not increase to the same extent. Therefore, iron will continue to 
penetrate or diffuse into the TiC porous media and alumina flow 
stops. The results in Figure 6 indicate that the alumina particles are 
expected to diffuse into TiC at a depth between 15 to 45 |im, depend- 
ing on particle size. The FESEM micrograph in Figure 1 (a) and its 
corresponding EDS elemental analysis (Table 1) strongly support the 
simulation MSPH results. It is thus clear why alumina accumulated 
at the innermost layer relative to the rotation axis, but the solid 
alumina layer isolated the iron particles. However, the iron particles 
remained in liquid phase and were able to stay in the molten medium 
for a longer period than alumina (Figure 8). This phenomenon 
allows the iron particles to wet and infiltrate the TiC porous media, 
and form a stronger composite structure. Figure 8 (b) displays how 
the iron particles surround TiC in zone (III). Therefore, the above- 
mentioned simulation has excellent potential in controlling particle 
diffusion. 




Time (s) 

Figure 5 | MSPH estimates of the particle velocity gradient (PVG) in the 
radial direction for three different alumina and iron particle sizes. 



SCIENTIFIC REPORTS | 4 : 3724 | DOI: 1 0.1 038/srep03724 



4 



1.6 



1.4 



1.2 



0.8 



0.6 



g 

£ 0.4 



0.2 



-♦-Alumina (h=2.5e-3; dp=1.0e-4) 
Alumina (h=3.5e-3; dp2=1.2e-4) 
-Alumina (h=1.5e-3; dp3=0.7e-4) 
-^Iron (h=2.5e-3; dp=1.0e-4) 
-^Iron (h=3.5e-3; dp2=1.2e-4) 
-^Iron (h=1.5e-3; dp3=0.7e-4) 




0 0.2 0.4 0.7 0.9 1.1 1.3 1.5 1.7 1.9 2.1 2.3 2.5 2.7 2.9 

Time (s) 

Figure 6 | Instantaneous change in radial direction versus time of alumina and iron particles during deposition. 



As explained earlier, the total mass and volume of the solution 
domain were approximated based in the particle method. As the 
number of particles approaches infinity, the approximation error 
moves toward zero. Essentially, both real and approximated solution 
domains are identical, thus, the exact analytical and MSPH solutions 
must be the same. However, because the simulation process included 
a limited number of particles, simulation error was present. 
Therefore, for a small number of particles the approximation error 
is big given that the particles are large. In the case of size dp2, shown 
in Figure 3 and Figure 5, the particle size was 120 um, which is the 
largest among the simulated cases. Thus, it is anticipated that the 
simulation error would be the largest. This helps clarify the signifi- 
cant difference in the two models' results. Moreover, the viscosity 
equation was also approximated based on particle size. Similarly, for 
particle size dp2, the viscosity approximation error was the largest. 
Especially at 0.9 and 2.1 seconds from the beginning of the reaction, 
the approximation error was quite significant (Figure 3). Given that 
at these two points abrupt changes in viscosity value occurred, par- 
ticle approximation is presumed to be even poorer. 

The MSPH simulation results are in agreement with the experi- 
mental observations. Alumina viscosity increases faster than iron, 
and alumina protects iron from external heat exchange, the iron 
remains in the molten medium longer than alumina. This phenom- 
enon allows the iron particles to wet and infiltrate the TiC porous 
media and form a strong composite. Mathematical modeling of par- 
ticle deposition velocity and viscosity can describe the process of Fe 
infiltration versus time. Melt viscosity increment due to the decrease 
in temperature leads to a significant speed reduction of particle 
deposition onto TiC. The MSPH method provides a very good 
estimation of velocity and velocity gradient of iron and alumina 
molten particles. The MSPH estimates the magnitude of molten iron 
and alumina infiltration, the depth, velocity and velocity gradient. 
The velocity gradient of particles reveals that the product velocity 
decelerates as the viscosity increases versus time. The prediction is in 
agreement with the experimental results, microstructure, and EDS 
observations. 

Methods 

Experimental methods. A centrifugal thermite (CT) machine facilitated rapid 
centrifugal acceleration as well as the temperature increase during the experiment. A 



bilayer graphite- steel compacted mold was fixed inside the CT reaction chamber. A 
high performance infrared thermometer, Raytek MM1MHSF3L, recorded the 
experimental real-time temperature data. The detailed procedure of using the 
centrifugal machine for thermite processing in metallic pipes is explained in recently 
published literature 31 . 

The starting materials, namely Al (<75 urn, 99% purity, Sigma Aldrich) and Fe 2 0 3 
(<5 urn, 97% purity, Sigma Aldrich) powders were dried for 8 hours and mixed in a 
Retsch PM 200 Planetary Ball Mill for 4 hours at 30-minute intervals. The green 
powder stoichiometric mixture was prepared according to the following reaction 
equation (l) 32 . 



Fe 2 0 3 +2Al^Al 2 0 3 +2Fe, AH= -836 kJmol~ 



(1) 



Titanium (Sigma Aldrich, — 100 mesh, 99.7% purity) and carbon (Sigma Aldrich, 
— 1000 mesh, 99.9% purity) elemental powders were dried for 7 hours at 125°C, 
following which they were mixed using a Planetary Ball Mill (Retsch PM 200) 
according to the reaction equation (2): 



Ti + C^TiC, AH = 



-183 kjmor 



(2) 



A 75 mm long carbon steel pipe with inner and outer diameters of 69 mm and 
75 mm, respectively, served as a bulk holder for graphite. The dried Ti + C green 
powder was compressed into a pellet. The pellet was then inserted into the graphite 
mold's pre-defined location as shown in Figure 7. The pellet was positioned near the 
head of the tube because of the higher thermal gradient 33 . 

The thermite mixture was fed into the tube. The rotation speed was subsequently 
increased to 280 g acceleration, followed by igniting an electric arc with a graphite 
electrode at the entrance of the tube. The onset of the thermite reaction was trailed by 
the titanium and carbon reaction, which occurred in response to the heat generated by 
the thermite (first) reaction. The released temperature was recorded by an infrared 
pyrometer. The detailed process and characterizations are explained extensively in a 
recently published paper by the authors 8 . 

Modeling and simulation of SHS using modified smoothed particle 
hydrodynamics (MSPH). Implementation of MSPH in the centrifugal thermite 
method. Figure 8 shows a schematic that represents the centrifugal SHS components 
in particle approximation. The TiC porous pellet, A1 2 0 3 and Fe Molten particles are 
shown just after the SHS reaction started (Figure 8 (a)). The representation of the 
infiltration of molten alumina and iron particles into the TiC pellet are seen in 
Figure 8 (b) and (c), which show the TiC pellet just after particle infiltration stopped 
and the particles approached the inner surface of the steel pipe. In this model, all 
particles are distinct in volume. 

Comprehensive MSPH formulations are provided in literature 26 ' 34 . The MSPH 
method was adapted to the TiC-Fe-Al 2 0 3 centrifugal SHS process. Cylindrical 
coordinates (r,0,z) were used to analyze the axisymmetric motion of particles, since 
the reaction occurred inside a cylindrical steel pipe. The function, f, describing a 
displacement variable, is independent of the angular position, 6, of a particle 
point. Particle approximation using MSPH is mainly dependent on the Taylor 
series expansion of function /(r,z) about point (r^), which is given as: 
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A: Before Reaction B: After Reaction 



Figure 7 | Schematic setup of the centrifugal thermite-assisted Ti-based functionally graded coating process; (A) before reaction, (B) after reaction. 



(3) 



where the derivatives are assessed at point x = (r^Zj) = (xi, yj. Neglecting the third 
and higher order derivative terms on the right side of equation (3) and multiplying 
both sides of the equation with a positive- value kernel function W(x — £, h) of 
compact support, or a smoothing length of 2 h, determines the number of particles 
involved in the approximation. The center of the particle position, x, its first and 
second derivatives, then integrate the resulting equations over the solution domain. 
The integration domain can be replaced by the compact support of the kernel func- 
tion since the kernel function is zero when the support is above 2 h. In this way, it is 
possible to estimate and write particle position, particle velocity, and particle velocity 
gradient in matrix form: 



E <p x Wijrrij/pj E (PiWijmj/pj 
j=i j=i 

N N 

E PiWijjcntj/pj E (Pi w ij, x m j/Pj 
j=i j=i 

N N 

E ViWijjmj/pj E (PiWijymj/pj 
j=i j=i 

N N 

E PiWijjcxmj/pj E (PiWij^mj/pj 
j=i j=i 

N N 

E (PiWij^rrij/pj E (PiWij^mj/pj 
j=i j=i 

N N 

E (PiWij^mj/pj E (PiWij^mj/pj 
j=i j=i 



N 

E? 

; = i 


eWijmj/pj 


N 

E? 

i=i 


eWij^mj/pj 


N 

E? 

i=i 


eWi^mj/pj 


N 

E? 


eWij^mj/pj 


N 

E? 

i=i 


eWy^rrij/pj 


N 

E? 

; = i 


eWij^mj/pj 



BF = LorB, l F,=L,I = 1,2,3, 



(4) 



Where 







jtfiWijmj/pj 
;=i 


( fi ) 




Zm^rrij/pj 
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■f 

J XXI 

fyyi 


>■ L=< 
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<fxyi . 




N 
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EfiWij^mj/pj 
w"=i ) 



and 



Pi = l, (p 2 =Xj-Xi, (p 3 = yj - yi , 



(p 5 = 



To obtain the solution for F, matrix B should be non- singular. This condition is 
already satisfied since the functions (p { are linearly independent and the number of 
particles in the compact support of the kernel function W for particles is at least 
6. All derivatives of kernel function W appearing in matrix B, must not be 
constants. Therefore, to easily approximate and calculate the radial positions and 
infiltration rates of iron and alumina into TiC in the centrifugal-assisted SHS 
method, the revised Gauss function W (x-^, h) in equations (4) is used as the 
kernel function. 
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stopped 



Figure 8 | Schematic of a TiC pellet and the molten alumina and iron particles at different positions. 
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No 



Solve the heat diffusion term 
in equ (6) by MSPH 



Calculate particle velo city 
and velocity gradient, equ (7) 



March the solution in time based on forward 
difference scheme, and update temperate and 
viscosity 




Figure 9 | MSPH computing method to predict particle properties. 



1.10081 



W(x-Z,h)=\ (hy/%) 2 
10 



-i\ z /h 2 



\x-£\<2h, 
\x-£\>2h. 



(5) 



In the present SHS model, smoothing length, or interaction radius (h), is selected 
such that on average 15-20 particles are interacting 26 . 

Solving the heat balance equation. Based on the geometric configuration of the SHS 
process, the temperature distribution is symmetric about the z-axis. Therefore, the 
thermal condition problem is simply a two-dimensional axisymmetric problem. For a 
unity volume, the heat balance equation or the transient heat flow governing equation 



Table 3 | Thermophysical properties of Fe-A^Oa-TiC particles 

Reaction 1 [equation (1 )] Reaction 2 [equation (2)] 





Fe 




Al 2 0 3 


TiC 


CJ kg" 1 K" 1 


473 




880 


881.98 


m, kg 


9.19E-02 




8.39E-02 


3.30E-03 


Po/ kg m" 3 


7800 40 




3950 41 


49 3 0 42 


p m/ kg rrr 3 


27 00 43 




2750 44 




K, w/m.K 


47 








C P/ kJ/kg.K 


0.49 








Q, kJ mol" 1 




145 




560 


q, kJ mol" 1 




836 




183 45 



dT 

"Tt" 



fd 2 T idT d 2 r 

V dr 2 r dr dz 2 



(6) 



where T is the transient temperature, t the time during the process, k is thermal 
conductivity, p is density, C p is the heat capacity, and q is the rate of heat generated 
during the reaction process. The last term in equation (6) functions as a heat sink from 
the reaction zone in the SHS process. The heat is consumed in preheating the reacted 
mixture. With the help of the initial conditions, the MSPH method can be utilized to 
solve equation (6). The starting mixture temperature and activation energy are given. 
Equation (6) is written for particle i and the values of T rri , T ri and T zzi , which are 
similar to T xxi , T xi and T yyi , and are determined by solving the linear simulation 
equations (4). 



fxi 



dT(r h Zi) _ 8T(xi,yi) 



BnL^J^BnW^mj/pj.Tj 



-- B M Li = ^ B ^ Wij, xx rrij / pj . 7) 
j=i 



-- B I5 L Y = B I5 Wij^nij/pj . Tj 



The solution is marching forward in time by the conditionally stable forward-dif- 
ference scheme. Upon obtaining the temperature distribution, it is possible to cal- 
culate the viscosity, particle velocity and velocity gradient. 

Solving the particle velocity and velocity gradient equations. Form Stoke's formula for 
particle motion in fluid, as well as normal gravity and equilibrium of forces on a 
particle in the SHS process under centrifugal acceleration 8 , the radial particle velocity 
V (T, t) as a function of temperature and time is given as: 



dr 


dx 


d 2 T(ri,Zi) 


d 2 T(x h yi) 


dr 2 


dx 2 


d 2 T{r uZi ) 


d 2 T(x hyi ) 


dz 2 


dy 2 



v(T C9 ty. 



.d 2 { Po - 



■p m )(39ARn 2 +g) 
18#7 



rj(T c ) = rj 0 exp 



Q_ 
R?-T c 



(7) 



(8) 



where R is the current radial position of the particles, rj (TJ is the current viscosity of 
the molten metal, Q is total activation energy, and the gas constant is R g . The rotation 
speed per second is n, g is the acceleration of gravity, and T c is the temperature. Particle 
size is d whereas p 0 and p m are the particle density and melt density, respectively. 

Since the determinant of matrix B is nonzero, its inverse, B~ l > exists. Thus, equation 
(4) can be rewritten as: 



F = B~ 1 T or F J = B~ 1 T I , J= 1,2,3, 



(9) 



Equation (8) gives the values of function/ (velocity, V) and its first-order and second- 
order derivatives at point x = (x i5 yi) in terms of velocity values at points neighboring 
x. Based on equation (9), MSPH approximation of radial velocity and radial velocity 
gradient can be explicitly written as: 



fi = V{r u zi) = V(x i ,y i )=B Il T l = J2" =l Bn W^/pj.Vj 



(10) 



f* 



dV{ri,Zi) _ dV(x h yi) 



il- 



ex 



-B I2 Ti = X^ =1 BnW^mj/pj.Vj (11) 



It is evident that for the MSPH method, kernel approximation of the velocity gradient 
is given in terms of the integration of the first-order derivative of W(x-cj, h) over its 
compact support. 

Coding and simulation procedure. MATLAB, Mathworks ® version 2012a, was 
employed to write the MSPH for solving equations (6) and (7). A flowchart of the 
solution procedure is shown in Figure 9. First, using the initial condition, the heat 
condition equation is solved and once the temperature distribution is known, vis- 
cosity, particle radial velocity, and radial velocity gradient are calculated. The solution 
marches forward in time by the conditionally stable forward-difference scheme. The 
time increment employed is 1 |is. Interaction radius or smoothing length (h) is 
chosen such that 15-20 particles are interacting on average 26 , where the total number 
of particles is 29. The MSPH simulation parameters are listed in Table 1. 

MSPH simulation of iron and alumina particles was conducted with parameters 
extracted from existing literature. Molten iron r\ 0 = 0.0065 kg/m.s 35 , total activation 
energy, Q = 706 kj/mol of Al-Fe 2 0 3 and Ti-C 3637 , gas constant R g = 8.31441 J/ 
K.mol 38 , and r\ Q = 0.052 N.s/m 2 is the liquid viscosity of alumina at 2408 K 39 . The 
other thermophysical properties are listed in Table 3. 
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